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JE 要 : 在 科尔沁 沙 地 采用 ”两 行 一 带 " 种 植 模式 低 密度 栽培 樟 子 松 进行 生态 修复 后 ,林地 土壤 水 分 动态 及 其 对 降 


雨 的 响应 影响 了 同类 型 地 区 是 否 可 以 持续 使 用 樟 子 松 进 行 
水 分 动态 特征 ,本 研究 综合 运用 原 位 观测 数值 模拟 的 方法 


雨 -土壤 水 分 响应 关系 。 结 果 表 明 :(1) 榜 子 松 信 工 林 显著 改变 了 地 区 水 分 分 布 , 裸 沙 地 2.0 m 处 深层 渗 漏 量 占 降 雨 


量 的 44.16% , 而 樟 子 松林 地 深层 渗 漏 量 仅 占 降雨 量 的 0.7%。 


生态 修复 。 为 了 研究 科尔沁 沙 地 南 缘 植被 修复 后 土壤 
,基于 土壤 水 分 实测 数据 校正 Hydrus- 1D 模型 ,探究 降 


(2) 监测 期 内 ,0.4 m 深 度 以 下 土壤 水 分 对 小 雨 无 响应 ， 


prz 


土壤 水 分 对 中 雨 的 响应 次 度 可 达 1.0m, 对 大 十 和 暴雨 的 响应 深度 涉及 整个 观测 剖面 。 随 着 土壤 深度 的 增加 ,水 分 


波动 幅度 呈现 减 小 的 趋势 。(3) 降雨 量 和 深度 较 浅 的 土壤 体积 含水 量 之 间 存 在 较 强 的 相关 关系 ,周期 为 周 .半月 的 
累计 降雨 量 与 各 层 土 壤 体积 含水 量 显著 相关 ;降雨 量 大 于 50 mm 时 ,能 保证 对 2.0 m 处 土壤 水 分 的 补给 .(4) 模型 的 
决定 系数 范围 在 0.61~0.85, 均 方 根 误差 范围 在 0.0061~0.0096 cm? * em? ,能 较 好 地 模拟 研究 区 土壤 水 分 的 动态 变化 


特征 , 且 深 层 模拟 精度 高 于 浅 层 。 人 研究 结果 对 科尔沁 沙 地 


意义 。 


雨 养 型 植被 造林 、 生 态 恢 复 和 水 资源 管理 具有 重要 的 


关键 词 : 降雨 人 渗 ; 土壤 水 分 ; 深层 渗 漏 ; 樟 子 松 ; 科尔沁 沙 地 


中 国 在 过 去 40 a 间 开展 了 广泛 的 植被 绿化 项 
目 ,森林 的 增加 大 幅度 增强 了 生物 地 球 化 学 和 生物 
物理 对 气候 的 影响 。 同 时 ,人 工 建 植 促进 植被 恢复 
可 显著 影响 土壤 化 学 计量 特征 ,也 对 土壤 物理 改良 
具有 积极 作用 ,有 效 改善 土壤 理化 性 状  。 在 干旱 
半 干 旱地 区 的 陆地 表层 系统 中 ,地 区 水 资源 平衡 受 
降雨 . 蒸 散 发 .径流 和 土壤 储 水 量 控制 ,其 中 土壤 水 
是 水 资源 重要 的 存在 形式 ,是 植被 恢复 格局 和 过 程 


为 全 国 沙漠 化 土地 最 集中 的 省 区 之 一 ,理解 植被 恢 
复 后 土壤 水 分 的 动态 变化 对 于 后 续 植 被 恢复 极其 
重要 。 

学 者 对 旱 区 人 工 植被 土壤 水 分 的 时 空 动 态 特 
征 及 其 再 分 配 开 展 了 大 量 研究 。Zhang 等 “指出 植 
被 恢复 后 ,干旱 荒漠 化 地 区 土壤 水 分 平衡 发 生 了 显 
著 变 化 。Cheng 等 “在 干旱 区 开展 了 大 量 关 于 人 工 
林 和 农田 的 深层 渗 漏 研究 ,发 现 人 类 活动 改变 了 沙 


的 主要 驱动 力 与 重要 生态 因子 。 降 雨 转换 而 来 
的 土壤 水 分 对 半 干 旱地 区 植被 恢复 与 重建 影响 显 
车 , 雨 养 型 植被 是 中 国防 沙 治 沙 工程 的 主要 植被 类 
型 ,其 高 度 依赖 土壤 水 生存 ”。 干 旱 少雨 水 资源 匮 
Z .植被 退化 以 及 沙漠化 加 剧 等 因素 导致 内 蒙古 成 
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量化 旱 区 水 资源 分 配 ,进行 生态 恢复 和 退耕 还 林 提 
供 了 科学 依据 。 目 前 在 旱 区 ,降雨 水 分 运 移 研 究 中 
一 般 通 过 对 相关 参数 进行 具体 的 修正 ,形成 特定 条 
件 下 的 水 分 运动 模型 学 。 其 中 ,土壤 水 分 模拟 是 目 
前 Hydrus-1D 应 用 最 为 广泛 的 领域 ,能 较 好 的 模拟 
干旱 半 干 旱 条 件 下 的 土壤 水 分 平衡 RIBUS ICI, 
研究 下 渗 对 降雨 格局 的 响应 。 

樟 子 松 是 一 种 适宜 于 中 国 北方 的 树种 , 节 水 耐 
盐 。 在 科尔沁 沙 地 ,降雨 是 生态 系统 和 土壤 水 分 最 
为 主要 的 水 分 补给 来 源 , 且 该 地 区 的 降雨 条 件 接近 
樟 子 松 引 种 区 的 降雨 下 限 , 因 此 樟 子 松 被 大 量 用 于 
营造 防风 固沙 林 来 治理 生态 环境 问题 "1。 在 生态 
脆弱 的 奈 曼 沙 区 种 植 樟 子 松 可 以 有 效 改善 土壤 结 
构 ,适宜 在 该 地 固沙 工作 中 推广 拉 。 采 用 ”两 行 一 
带 ” 的 种 植 模式 低 密度 栽培 樟 子 松 进行 生态 修复 
后 ,林地 土壤 水 分 动态 及 其 对 降雨 的 响应 研究 极其 
重要 ,这 影响 了 同类 型 地 区 是 否 可 以 持续 使 用 樟 子 
松 进行 生态 修复 。 但 现 阶 段 对 科尔沁 地 区 “两 行 一 
带 ” 式 樟 子 松 的 土壤 水 分 平衡 研究 进展 不 大 ,因为 
深层 土壤 水 分 运 移 过 程 不 清晰 导致 没有 完整 研究 
样 地 尺度 的 地 上 地 下 部 分 水 量 平衡 。 因 此 ,本 研究 
以 樟 子 松 为 对 象 ,采用 野外 控制 定位 试验 监测 法 和 
数据 模拟 分 析 法 相 结合 ,分析 土壤 水 分 时 空 动 态 及 
不 同 深 度 土壤 水 分 对 降雨 的 响应 ,论证 Hydrus-1D 
模型 对 研究 区 土壤 水 分 运动 模拟 的 可 靠 性 ,探究 使 
用 樟 子 松 进 行 植被 恢复 后 地 区 土壤 水 分 的 动态 
特征 。 


1 数据 与 方法 


1.1 研究 区 概况 

人 研究 区 位 于 通辽 市 奈 曼 旗 (120 19' 40”~ 
121^35'40" E, 42^ 14'40"-43* 3220" N) ZR JE E A Zi pf 
林场 (图 1) ,科尔沁 沙 地 南 缘 。 该 地 区 海拔 为 240~ 
700m ,地势 西南 高 ,东北 低 "”。 属 于 温带 大 陆 性 季 
风气 候 ,年 平均 气温 7.0"C ,年 平均 降雨 量 354.3 
mm ,降雨 主要 分 布 在 5 一 10 月 。 夏 季 高 温 降雨 较 
多 ,冬季 低温 且 降 雪 较 少 ,昼夜 温差 大 , 日照 时 间 
长 。 该 地 区 风沙 土 广 布 ,土壤 孔 际 大 ,透水 性 较 
强 。 同 时 还 分 布 有 草 铅 土 母 质 , 蓄 水 保 肥 能 力 强 ， 
通气 性 较 差 ,不 耐 旱 涝 中 。 地 带 性 植被 为 温带 草原 
带 , 次 生性 明显 。 以 半 干 旱地 区 的 植物 种 为 主 , 天 
人 然 分 布 的 植物 种 主要 有 山 奉 (Armeniaca sibirica), H 
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图 1 研究 区 地 理 位 置 
Fig. 1 Location ofthe study area 


F (Imperata cylindrica ) KEFF (Stipa grandis) FE 
(Leymus chinensis ) H (Salix gordejevii ) 等 |。 

L2 野外 调查 与 数据 采集 

1.2.1 试验 地 选择 ”选择 奈 林 林场 地 势 平缓 的 榜 子 
PA (Pinus sylvestris var. mongolica ,PSM) 林 地 ,以 减少 
坡度 对 研究 结果 的 影响 。 研 究 区 樟 子 松 已 种 植 32 
a, 采 用 行 带 式 种 植 ,两 行 之 间距 离 20 m。 样 地 为 长 
方形 ,大 小 300 mx500 m, 内 有 2500 ~ 3500 株 樟 子 
松 ,坡度 缓 , 视 为 平地 。 土 壤 类 型 为 沙土 ,无 地 表 径 
流 , 地 下 水 平均 埋 深 5.75 m。 在 行 带 完整 的 样 地 里 ， 
选取 胸径 中 位 数 植株 作为 试验 对 象 , 对 样 地 200 株 
榜 子 松 进行 调查 后 各 项 平均 指标 如 表 1 所 示 。 为 了 
避免 空间 异 质 性 带 来 的 试验 误差 ,选取 条 件 相似 的 
3 块 试 验 标准 区 域 作为 重复 试验 样 地 ,避免 因 植物 
种 群 不 同 带 来 的 试验 差异 ,3 块 样 地 之 间距 离 大 约 
100 m。 同 时 选取 一 个 相 邻 的 裸 沙 地 作为 对 照 组 ， 
先 将 周围 树 根 挖 断 ,再 使 用 0.5 mm 塑料 膜 隔离 出 一 
个 周 长 20 m, T4 m 的 独立 土壤 空间 ,防止 周边 根系 
侵入 影响 降雨 水 分 和 人 渗 。 试 验 场地 设计 如 图 2 
所 示 。 

1.2.2 数据 采集 方法 ”人 研究 基于 原 位 试验 ,需要 获 
取 降 雨量 .土壤 含水 量 ,深层 渗 漏 量 等 信息 ,需要 获 
取 的 参数 如 下 : 


表 1 试验 样 地 樟 子 松 各 项 指标 平均 值 
Tab.1 Average value of each index of PSM 
in the study area 
树种 树 高 /m 枝 下 高 lm。 ”胸径 /cm ”胸围 /cm? 
樟 子 松 10.33 7.63 7.22 160.82 
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4» 4» (数据 记录 器 ) d> 4» Kh) JgAETRURISEZK lem d) ; h 2g FRZJZK3k Cem) ; 
piss "ER Sh) 为 根系 吸水 速率 (em.d-))。 


毛管 水 平衡 部 分 深层 渗 漏 仪 深层 渗 


裸 沙 地 图例 试验 标准 地 
一 土壤 水 分 传感器 
| 降雨 水 分 
图 2 试验 地 设计 


Fig. 2 Experimental Site design 


(1) 降雨 量 : 雨 量 计 的 安装 参照 标准 气象 站 安 
装 规范 ,在 样 地 安装 两 套 地 表 雨 量 监测 系统 ,统计 
降雨 的 雨量 .次数 .历时 等 指标 ,监测 数据 每 10 min 
记录 一 次 。 采 用 美国 AVALON 公司 生产 的 AV- 
3665R 型 雨量 传 感 吕 ,分辨 率 0.2 mm ,精度 +2%。 

(2) 土壤 水 分 :采用 土壤 水 分 传感器 分 层 监测 
样 地 2 m 内 土壤 的 水 分 变化 。 根 据 现场 挖掘 的 土壤 
剖面 ,结合 试验 地 草本 根系 特征 和 土壤 分 层 情况 ， 
将 樟 子 松林 地 土壤 划分 为 8 个 土壤 水 分 观测 层 ,并 
在 0.4m.0.6m.0.8m.10m.1.2m.1.4m、1.7m 和 
2.0m 深 度 处 安装 土壤 水 分 传感器 ,将 裸 沙 地 土壤 划 
分 为 6 个 土壤 水 分 观测 层 , 并 在 0.3 m 0.6 m,0.9 m, 
1.2 m,1.5 m 和 2.0 m 深 度 处 安装 土壤 水 分 传 感 絮 ， 
对 比 造林 后 降雨 水 分 再 分 配 过 程 受 到 植被 影响 而 
产生 的 差异 。 监 测 数据 每 10 min 记录 一 次 。 

(3 ) 深层 渗 漏 :使 用 深层 渗 漏 仪 对 样 地 2.0 mm 处 
深层 渗 漏 进行 连续 监测 ,该 仪器 利用 毛管 水 力 提升 
高 度 设 计 一 个 平衡 部 , 待 测 土屋 有 水 分 流量 进入 平 
衡 部 后 在 仪器 下 部 产生 渗 漏 ,从 而 直接 测量 深层 渗 
漏 量 。 在 安装 渗 漏 计量 仪 时 为 避免 扰动 观测 层 土 
壤 , 采 取 侧 挖 埋设 ,减少 对 原 位 土壤 层 的 干扰 ,保持 
原 位 土壤 的 完整 性 。 仪 器 计量 分 辩 率 为 0.2 mm, 精 
度 +2%。 

1.3 模型 计算 
1.3.1 水 流 方程 ”采用 Hydrus-1D 模 型 对 樟 子 松林 
地 土壤 水 分 运 移 过 程 进行 模拟 ,水 分 以 垂 向 运动 为 
E RHEE ELA IB Richards 7i fe ^ f$ : 

20-2 KOE $ 中 -s(h) (1) 
WP: 0 为 土壤 体积 含水 量 (cm - em?) ; 为 时 间 
(d); z 为 垂 向 坐标 (em) ,假设 地 表 为 0, 向 下 为 正 ; 


土壤 水 分 特征 曲线 和 非 饱和 导 水 率 用 van Ge- 


nuchten-Mualem 方程 拟 合 : 


ch ^ oun 
60-1 — [ren] (2) 
0 ， hz 0 
A^ 
| " 
K(h)- ZI , h«0, m=1-1/n (3) 
K, h>0 
0-0 
s-( ) (4) 
(6.- €.) 


式 中 : 0, 9 E BESEAR AKE Cem em); 0, 为 土壤 饱 
和 含水 量 (cm em?) ; KK 为 饱和 导 水 率 (cm d"); S, 
为 无 量 纲 的 相对 含水 量 ; a um n LOU C 
1.3.0 根系 吸水 ”Hydrus-1D 将 根系 吸水 转化 为 水 
分 胁迫 函数 ,提供 Feddes 和 S-shape 两 种 根系 吸水 
模型 。 本 研究 采用 Feddes 模型 来 模拟 模子 松林 地 
的 根系 吸水 过 程 ,具体 表达 式 如 下 : 

S(h) 2 a(h)S, 2 o(h)b(x)T,, (5) 
式 中 : SQ) 为 根系 吸水 速率 (cmd '); alh) 为 水 分 胁 
迫 反 应 系数 ; 5, 为 潜在 吸水 速率 (cmd'); Ha) 为 根 
系 吸水 分 配 密度 函数 ; T, 为 潜在 蒸腾 速率 (cmd')。 
1.3.3 时 间 和 空间 离散 ”由 于 降雨 集中 在 7 一 9 月 
且 此 时 段 为 樟 子 松 生 长 旺盛 的 季节 ,因此 选择 2021 
年 6 月 1 日 一 10 月 31 日 为 模拟 时 段 ,时 间 离 散 单元 
为 d, 共 153 d。 采 用 时 间 离 散 化 的 方式 进行 模拟 ， 
设置 起 始 时 间 为 0, 模 拟 结 束 时 间 为 153 d, 初 始 时 
间 步 长 为 0.001 qd, 最 小 时 间 步 长 为 0.001 d, 最 大 时 
间 步 长 为 1d。 模 拟 土 层 厚 度 为 2m, 以 0.01 m 离 散 
单元 将 其 分 割 为 200 个 单元 ,垂直 离散 为 201 Ai 
点 。 根 据 土壤 性 质 ,将 土壤 剖面 模型 分 为 8 层 ,设置 
8 个 观测 点 。 
1.3.4 边界 及 初始 条 件 ”模型 上 边界 条 件 设置 为 可 
积 水 的 大 气 边界 ,下 边界 设置 为 自由 排水 边界 。 洪 
Tr H z& PS OH] FAO 推荐 的 Penman-Monteith 公式 
计算 ,Hydrus-1D 自动 处 理 蒸 腾 量 每 日 变化 。 
1.3.5 模型 参数 率 定 与 评价 指标 ”模型 水 力 参 数 的 
准确 性 对 一 维 土 壤 水 分 运 移 模 型 能 否 反 映 实际 情 
况 具有 重要 的 决定 性 ,因此 需 对 参数 9, 、0,、a un. 
K, 进行 不 断 地 优化 ,使 其 更 加 符合 实际 情况 。 参 数 
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优化 时 Hydrus-1D 使 用 结合 了 牛顿 法 和 最 速 下 降 法 
的 LM 法 ,进行 目标 函数 的 最 小 化 ,并 为 优化 参数 提 
供 置 信 区 间 '。 本 研究 通过 不 同 降雨 条 件 下 的 实 
测 土壤 体积 含水 量 对 参数 进行 校正 与 优化 ,获取 的 
参数 见 表 2。 

选取 决定 系数 ( 尺 ) .相对 误差 (RE) 、 均 方 根 误 
差 (RMSE) 和 Nash 效 率 系数 (NSE) 作 为 评价 指标 进 
TT 25 41 Br ^" ,衡量 模型 的 模拟 效果 与 评价 模型 精 
度 。 其 中 ,RE 反映 模型 模拟 过 程 中 模拟 值 与 实测 值 
之 间 的 相对 误差 ;RMSE 反映 模拟 值 与 实测 值 的 绝 


对 误差 的 平均 程度 , 值 越 接 近 于 0, 模 拟 精度 越 高 ; 


NSE 表示 模拟 效果 随时 间 变 化 的 符合 程度 , 值 越 接 
近 1 模 拟 效 果 越 好 ”1。 
1.4 数据 处 理 

统计 发 现 2021 年 4 一 10 月 降雨 量 约 占 全 年 降 
雨量 的 98%, 故 选择 2021 年 4 月 1 日 一 10 月 31 日 作 
为 监测 期 分 析 降 雨水 分 分 配 。 将 24h 降 雨量 分 为 4 
个 等 级 :小 雨 (0.1~9.9 mm) , "PR (10.0-24.9 mm) , 


X Bj (25.0-49.9 mm) , a& M (50.0-99.9 mm) ™® ~ fii 


用 水 量 平 衡 方程 及 人 渗 系 数 公式 ,初步 计算 研究 区 
樟 子 松林 地 、 裸 沙 地 降雨 水 分 的 分 配 ,具体 表达 式 


如 下 : 


P-D=ASWC+ET (6) 
-D 
R=5 (7) 


式 中 : P HENE (Imm) ; D 为 深层 渗 漏 量 (mm ) ; 
ASWC 为 土壤 含水 量 的 变化 值 (mm); ET 2928 IU 
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量 (mm); R 为 人 渗 系 数 。 


应 用 Excel、 Matlab R2021b 软件 对 数据 进行 整 


理 分 析 , 并 采用 IBM SPSS Statistics 26 软件 进行 
Pearson 相关 性 分 析 , Origin 2022 软件 绘图 ,Hydrus- 


ID 建 立 模 型 。 


2 结果 与 分 析 
2.1 林地 与 裸 沙 地 降雨 水 分 分 配 


从 表 3 可 知 , 樟 子 松林 地 降雨 量 共 572.6 mm, T 
层 渗 漏 系数 为 0.70%。 其 中 ,0~2.0 m 土壤 中 储存 降 
雨水 分 37.34 mm, fj 4.0 mm 的 降雨 水 分 渗 漏 到 2.0 
m 以 下 ,其余 531.26 mm 水 分 为 植被 截留 RAE 
量 。 虽然 , 裸 沙 地 2 m 人 处 深层 渗 漏 量 大 于 樟 子 松林 
地 深层 渗 漏 量 , 但 林地 植物 根系 利用 水 分 且 蒸 散 发 
量 损失 较 大 ,监测 期 裸 沙 地 蓄 水 量 明显 多 于 林地 蓄 
水 量 。 使 用 模子 松 进行 生态 修复 后 改变 了 降雨 水 
分 的 分 布 , 但 樟 子 松林 地 仍 可 维持 大 气 降水 .地表 
水 .土壤 水 和 地 下 水 4 水 转换 ,并 保证 降雨 水 分 对 深 
层 渗 漏 或 者 地 下 水 的 补给 。 林 内 土壤 体积 含水 量 
随 着 降雨 量 的 增 减 而 呈现 相同 的 变化 趋势 (图 3) ,4 
月 .5 月 降雨 量 较 少 ,各 层 土壤 体积 含水 量 缓慢 增 
加 ;6 一 8 月 各 层 土壤 体积 含水 量 随 降雨 量 增加 而 迅 
速 上 升 ;9 月 降雨 量 减少 ,各 层 土壤 体积 含水 量 也 随 
之 降低 。 其 中 ,林地 降雨 主要 集中 在 7 一 9 月 ,但 2m 
以 下 深层 渗 漏 主要 发 生 在 5 一 7 月 ,8 月 深层 渗 漏 量 
仅 有 0.4 mm。 


表 2 研究 区 土壤 水 力学 参数 


Tab. 2 The soil hydraulic parameters in the study area 


TEE EE /m 0, K(cm^* cm?) 0,/(cm* em?) a / (cm!) n K./(cm*d^) l 
0~0.4 0.0178 0.415 0.0295 2.150 645.00 0.5 
0.4~0.6 0.0185 0.208 0.0310 2.510 815.00 0.5 
0.6~0.8 0.0017 0.305 0.0270 2.665 650.48 0.5 
0.8~1 0.0010 0.267 0.0217 2.780 736.00 0.5 
1~1.2 0.0010 0.271 0.0290 2.643 806.47 0.5 
1.2~1.4 0.0373 0.399 0.0300 2.580 852.00 0:5 
1.441.7 0.0246 0.398 0.0282 2.690 766.00 0.5 
1.7 ~ -2.0 0.0179 0.433 0.0187 4.630 1060.00 0.5 

注 : 0, 表示 残余 含水 量 ; 0, 表示 饱和 含水 量 ; K, 表示 饱和 导 水 率 ; a /表示 模型 参数 。 
表 3 试验 样 地 降雨 水 分 分 配对 比 
Tab.3 Comparison of rainfall distribution in experimental plots 
监测 期 样 地 类 型 ” 降雨 量 /mm ， 蓄 水 变化 量 /am ”2.0 m 以 下 渗 漏 量 /mm ” 燕 散发 量 /mm ” 深层 渗 漏 系数 /% 
20214E4—10 H 樟 子 松 572.6 37.34 4.0 531.26 0.70 
2021 年 4 一 9 月 裸 沙 地 571.6 88.70 2524 230.50 44.16 
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图 3 月 尺度 林内 降雨 量 ,深层 渗 漏 量 和 土壤 含水 量 动态 


Fig. 3 Dynamics ofrainfall, deep Soil recharge and soil 


water content at monthly scale 


2.2 林地 降雨 量 与 土壤 体积 含水 量 关 系 
2.2.1 林地 土壤 水 分 对 降雨 的 动态 响应 ”降雨 后 部 
分 雨水 下 渗 进 入 土壤 ,使 得 土壤 体积 含水 量 随 之 变 
化 。 从 图 4 可 以 看 出 ,在 雨季 不 同 深度 土壤 体积 含 
水 量 波动 频繁 ,每 次 波动 都 与 降雨 事件 相对 应 。 监 
测 期 共 记 录 77 次 降雨 事件 ,引起 林地 土壤 体积 含水 
量 明显 波动 的 共有 15 次 ,其 中 最 小 的 为 7 月 1 H 8.4 
mm 的 降雨 ,最 大 的 为 8 月 20 H 81.6 mm 的 降雨 。 不 
同 降雨 事件 引起 的 土壤 水 分 波动 情况 不 同 , 例 如 7 
H 9 H 45.8 mm 的 降雨 引起 了 8 个 观测 层 含水 量 的 
变化 ,对 0~2.0 m 土 壤 水 分 补给 作用 明显 ;而 9 月 6 
H 13.4 mm 的 降雨 只 引起 了 0~0.6 m 土 壤 水 分 的 波 
动 ,对 樟 子 松林 地 0.6 m 以 下 土 层 无 直接 补给 作用 。 
观测 层 中 ,林地 1.2 m 处 土壤 体积 含水 量 最 


0.4m 
-一 1.2m 
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图 4 樟 子 松林 地 不 同 深度 土壤 水 分 对 大 气 降雨 的 响应 
Fig. 4 Response of soil moisture at different depths to 
rainfall in PSM forest land 


定 ,4 月 林地 1.2 m 深 度 土壤 体积 含水 量 稍 大 于 1.4m 
处 土壤 体积 含水 量 ,在 5 一 7 月 两 处 土壤 体积 含水 量 
近乎 相等 。7 月 雨季 开始 ,1.2 mm 深度 土壤 体积 含水 
量 大 于 1.4 m 深 度 土壤 体积 含水 量 , 这 两 处 土壤 体 
积 含水 量 在 全 年 均 高 于 其 他 观测 点 的 土壤 含水 量 。 

由 樟 子 松林 地 各 观测 层 土 壤 体 积 含水 量变 化 
可 知 ( 表 4) ,8 个 观测 层 土壤 体积 含水 量 最 低 值 均 出 
现在 4 月 ,最 高 值 除 1.7 m 处 在 9 月 外 ,其 余 7 层 最 高 
值 均 出 现在 8 月 。 观 测 剖 面 4 一 10 月 最 大 土壤 体积 
含水 量 为 15.51%, 出 现在 8 月 20 日 1.2 m EEA, 
4 一 10 月 最 低 土壤 体积 含水 量 为 0.66% ,出 现在 4 月 
1 H 1 m-EJZAb, 0.4 m,0.6m,0.8m.,1.0m,1.2 m, 
1.4m、1.7m 和 2.0 m 处 土壤 体积 含水 量变 化 幅度 分 
别 14.03% 、13.00% 、13.64% 、10.13% , 9.8196 、 
9.3896 .8.66% 和 7.68%。 由 此 可 知 , 随 着 土壤 深度 的 
增加 ,水 分 波动 幅度 呈现 减 小 的 趋势 。0.4 m 深 度 土 
壤 体 积 含水 量 波动 最 大 ,受降 雨 影响 剧烈 。 
2.2.2 相关 分 析 ”大气 降雨 是 樟 子 松林 地 土壤 水 分 
的 主要 来 源 ,因此 ,降雨 量 与 不 同 深度 的 土壤 体积 
含水 量 之 间 可 能 存在 相关 关系 。 考 虑 到 降雨 后 渗 
漏 的 延迟 性 ,分别 统计 观测 周期 为 日 周 .半月 月 
的 累计 降雨 量 和 各 层 土壤 体积 含水 量 ,进行 Pearson 
相关 分 析 ,结果 如 表 5 所 示 。 对 于 0.4 m、0.6 m、0.8 
ml10m 和 1.2m 土 层 ,其 土壤 体积 含水 量 与 累计 降 
雨量 之 间 在 日 周 .半月 .月 尺度 上 均 存 在 相关 关系 
(P<0.05), 其 中 0.4m.0.6m.0.8m.1.0m 和 1.2m 土 
层 ,其 土壤 体积 含水 量 与 累计 降雨 量 在 日 、 周 .半月 
尺度 上 显著 相关 (P<0.01); 周 累计 降雨 量 与 1.4 m、 
1.7 m 土 层 体 积 含水 量 在 0.05 水 平 上 显著 相关 ,与 其 
余 土 层 体积 含水 量 在 0.01 水平 上 显著 相关 ;半月 累 
计 降 雨量 和 1.7 m 土 层 体积 含水 量 在 0.05 KF E m 
著 相 关 , 和 其 余 土 层 体积 含水 量 在 0.01 水 平 上 显著 
相关 。 

总 体 来 说 ,降雨 量 和 深度 较 浅 的 土壤 体积 含水 
量 之 间 存 在 较 强 的 相关 关系 。 周 期 为 周 .半月 的 累 
计 降 雨量 与 各 层 土 壤 体积 含水 量 显著 相关 (P< 
0.05) , 当 周 期 为 月 时 其 相关 性 整体 减弱 ,而 周期 为 
日 时 其 相关 性 只 在 较 浅 的 土 层 中 体现 。 
2.2.3 降雨 量 与 下 渗 深 度 分 析 由 试验 结果 可 知 下 
渗 深 度 到 达 0.4 m 处 的 最 低 降 雨量 为 13.4 mmo 3€ 
择 试 验 期 间 10 场 典型 独立 降雨 与 下 渗 深 度 进行 指 
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RA 樟 子 松林 地 不 同 深度 土壤 体积 含水 量变 化 峰值 
Tab. 4 Peak value of soil water content at different depths of PSM forest land 


2021 年 4 一 10 月 


最 小 值 /% 出 现 日 期 /月 -日 最 大 值 /% 出 现 日 期 /月 -日 差 值 /% 
0.4 1.29 04-01 15.32 08-20 14.03 
0.6 1.00 04-19 14.00 08-20 13.00 
0.8 0.67 04-01 14.31 08-20 13.64 
1.0 0.66 04-01 10.79 08-20 10.13 
1:2 5.70 04-01 15.51 08-20 9.81 
1.4 5.01 04-01 14.39 08-21 9.38 
1:7 1.90 04-01 10.56 09-21 8.66 
2.0 2.00 04-19 9.68 08-21 7.68 
表 5 降雨 量 与 不 同 深度 土壤 含水 量 的 相关 关系 
Tab.5 The correlativity between rainfall and soil water content at different depths 
时 间 尺 度 
0.4 m 0.6 m 0.8 m 1.0 m 1.2 m 1.4 m 1.7 m 2.0 m 
H 0.473” 0.502” 0.494” 0.410” 0.296" 0.173 0.100 0.124 
周 0.609™ 0.600” 0.595” 0.541” 0.622” 0.378 0.424 0.5347 
半月 0.735” 0.729” 0.724” 0.687” 0.775” 0.857” 0.578 0.684" 
0.811° 0.817 0.833* 0.824 0.876" 0.924" 0.68 0.870 


注 :* 表 示 在 0.05 级 别 ( 双 尾 ) 相 关 性 显著 ,** 表 示 在 0.01 级 别 ( 双 尾 ) 相 关 性 显著 。 


数 .S 型 方程 多项式 拟 合 分 析 , 最 优 拟 合 结果 如 图 5 
所 示 。 由 图 5 可知, 在 一 定 范围 内 ,下 渗 深 度 随 降雨 
量 的 增 大 而 增加 , 当 降 雨量 大 于 50 mm 时 ,能 保证 对 
2.0 m 处 土壤 水 分 的 补给 。 此 外 ,两 次 中 雨 雨量 分 别 
为 16.8 mm(6 月 19 日 ) 和 18.8 mm(8 月 24 日 ) ,数值 
接近 ,但 入 渗 深 度 分 别 为 0.6 m 和 1.0 m, 42 0.4 
m。 由 此 可 见 , 降 雨量 一 定时 ,下 渗 深 度 仍 可 能 存在 
较 大 差异 ,独立 降雨 的 入 渗 深度 还 受到 降雨 强度 、 
土壤 初始 含水 率 等 的 影响 。 
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图 5 降雨 量 与 下 渗 深 度 拟 合 
Fig.5 Fitting of rainfall and infiltration depth 


2.3 林地 土壤 水 分 运 移 模 拟 

在 参数 优化 后 ,根据 所 建 模型 模拟 得 到 土壤 水 
分 运 移 结 果 。 由 图 6 可 知 ,Hydrus-1D 能 够 模拟 降 
雨 信 渗 过 程 中 土壤 水 分 的 动态 变化 ,人 研究 区 土壤 水 
分 运 移 模 拟 值 和 实测 值 具有 基本 一 致 的 走向 变 
化 。 总 体 来 看 ,没有 降雨 发 生 的 6 月 和 10 月 拟 合 度 
较 有 降雨 时 的 7 一 9 月 拟 合 度 低 ,0.4 m、1.7 m 深 度 模 
拟 效 果 与 其 他 层 相 比较 差 , 其余 深 度 土壤 水 分 实测 
值 均 坐 落 于 模拟 曲线 附近 。 在 9 月 21 日 最 后 一 场 
降雨 后 ,0.4 m、1.4 m,2.0 m 处 模拟 结果 较 好 ,其 余 深 
度 土壤 水 分 被 低估 。 

使 用 实测 土壤 体积 含水 量 与 模型 模拟 值 验 证 
Hydrus-1D 在 研究 区 的 适用 性 ,检验 结果 显示 模型 
整体 决定 系数 为 0.876, 除 了 0.4m 处 尼 为 0.615 外 ， 
其 余 土 层 的 让 在 0.63~0.85 范围 内 , 拟 合 结果 较 好 。 
选择 具有 代表 性 的 4 个 观测 层 (0.4 m、0.8 m. 1.2 m 
和 2.0 m) 进行 分 析 , 下 文中 NSEis、 REG, R^; 与 
RMSE,; 分 别 表示 1.2 m 处 观测 层 的 Nash 效率 系数 、 
相对 误差 ,决定 系数 、 均 方 根 误差 ,其余 3 个 观测 层 
表示 方法 与 此 相同 。 各 土 层 土壤 水 分 模拟 效果 和 
评价 结果 如 图 7 所 示 ,4 个 观测 层 Nash 效率 系数 均 
大 于 0.5, 在 可 接受 范围 内 ;相对 误差 绝对 值 均 小 于 
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图 6 不 同 深度 土壤 含水 率 观测 值 与 模拟 值 比较 


Fig.6 Comparison of measured and simulated values of soil water content at different soil depths 


5%, 均 方 根 误 差 均 小 于 0.008, 模 拟 精 度 较 高 。 其 
中 ,0.8 m、1.2 m,2.0 m 处 观测 层 土壤 水 分 模拟 效果 
较 好 。 从 各 土 层 的 模拟 结果 对 比 来 看 ,Nash 效率 系 
数 表现 为 :NSE,o>NSE1;>NSEos>NSEw, 相 对 误差 的 
绝对 值 表 现 为 :RE<REus<REo<RE,。, 均 方 根 误差 
表现 为 :RMSE,uo<RMSE,,<RMSEu<RMSFEus ,决定 系 
BUCASHBION ROS RS Rss Ross 

综 上 可 知 , 经 参数 优化 后 的 Hydrus-1D 模 型 对 
研究 区 土壤 水 分 运动 模拟 有 较 好 的 适用 性 ,并 且 深 
层 土壤 含水 率 的 模拟 精度 高 于 浅 层 土壤 。 模 型 可 
通过 预测 的 降雨 来 模拟 土壤 水 分 的 变化 ,为 科尔沁 
沙 地 植树 造林 提供 科学 依据 。 


3 讨论 


3.1 林地 与 裸 沙 地 深层 渗 漏 分 析 

裸 沙 地 2.0 m 深 度 处 渗 漏 量 占 降雨 量 的 44.16%， 
樟 子 松林 地 渗 漏 量 占 降雨 量 的 0.7% ,使 用 樟 子 松 进 
行 生 态 修 复明 显 改变 了 降雨 水 分 的 分 布 。 深 层 土 
壤 水 的 补给 和 更 新 通常 比较 缓慢 ,在 干旱 半 干 旱地 
区 更 为 明显 ,因此 ,植树 造林 容易 造成 深层 土壤 干 
燥 化 现象 ,进而 阻止 上 层 土壤 与 地 下 水 之 间 的 水 交 
换 ,导致 深层 渗 漏 量 减少 下 。 林 木 会 对 降雨 转化 而 
来 的 土壤 水 分 进行 吸收 利用 , 旦 榜 子 松 作为 深 根 性 
植被 能 够 吸收 利用 深层 土壤 水 ,进一步 导致 降雨 补 
给 的 深层 渗 漏 量 减少 。 在 降雨 过 程 中 ,降雨 特性 也 
会 影响 植物 的 穿 透 和 截留 ,从 而 影响 到 降雨 水 分 的 


分 配 。Zhang 等 "指出 裸 沙 区 与 植被 恢复 区 土壤 芝 
发 量 差 异 的 主要 原因 是 生物 结 皮 的 形成 ,植被 恢复 
后 ,土壤 蒸 散发 随 之 改变 ,进而 影响 水 分 深层 渗 漏 
ti ,导致 植被 修复 后 形成 的 樟 子 松林 地 与 裸 沙 地 的 
深层 渗 漏 量 产生 较 大 的 差异 。 我 们 的 研究 表明 ,种 
植 模子 松 后 地 区 水 分 分 布 发 生 了 改变 ,因为 模子 松 
生存 耗 水 ,将 降雨 水 分 拦截 在 浅 层 土壤 中 并 消耗 ， 
减少 了 地 区 深层 渗 漏 。 虽 然 研 究 区 造林 后 降雨 水 
分 分 配 发 生 了 较 大 变化 ,但 2021 年 樟 子 松 林地 仍 可 
维持 4 水 转换 ,并 保证 降雨 水 分 对 深层 渗 漏 或 者 地 
下 水 的 补给 。 

此 外 ,监测 期 降雨 主要 集中 在 7 一 9 月 ,但 2.0m 
处 深层 渗 漏 主要 发 生 在 5 一 7 月 ,8 月 深层 渗 漏 量 仅 
有 0.4mm。 发 生 此 现象 是 因为 7 一 9 月 为 樟 子 松 生 
长 季 ,土壤 中 的 雨水 大 多 被 吸收 利用 , 且 此 时 的 闵 
散发 量 较 大 ,水 分 很 少 渗 漏 到 2 mm 以 下 。 而 4 一 5 月 
是 冻 融 季节 ,地 表土 壤 蒸 发 较 少 ,并 且 深 层 土 壤 和 
浅 层 土壤 存在 温度 差 ,导致 有 水 分 渗 漏 ,此 时 樟 子 
松 生长 活动 较 弱 ,吸收 利用 的 水 分 也 比 生长 季 少 。 
因此 ,在 7 一 9 月 樟 子 松 生长 旺 成 的 季节 ,即使 降雨 
量 增 加 ,水 分 大 多 被 林木 利用 , 深层 渗 漏 量 并 未 
增加 。 
3.2 林地 土壤 水 分 对 降雨 的 动态 响应 

樟 子 松林 地 0.4 m,0.6 m,0.8 m,1.0m,1.2 m, 
1.4m、1.7m 和 2.0 m 人 处 土壤 体积 含水 量变 化 幅度 分 
别 为 14.0396 , 13.0096 、13.64% 、10.13% 、9.81% 、 
9.3896 、8.66% 和 7.68%。0.4 m 深度 土壤 体积 含水 量 
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图 7 不 同 深度 土壤 水 分 运 移 模 拟 结果 验证 


Fig.7 Verification of simulation results of soil water transport at different soil depths 


波动 最 大 ,受降 雨 影响 剧烈 。 随 着 土壤 深度 的 增 
加 ,水 分 波动 幅度 呈现 减 小 的 趋势 。 降 雨 过 程 的 持 
续 使 士 壤 含水 率 逐 渐 增加 ,水 分 主要 在 重力 作用 下 
运动 ,深层 土壤 受降 雨 影响 小 ,水 分 逐渐 达到 稳 渗 
状态 ,因此 ,深层 土壤 含水 量 随时 间 的 持续 变化 比 
浅 层 小 。 将 24h 降 雨量 分 等 级 可 得 ,监测 期 内 0.4 m 
深度 以 下 土壤 水 分 对 小 雨 无 响应 ,土壤 水 分 对 中 雨 
的 响应 深度 可 达 1.0 m, 对 大 雨 和 暴雨 的 响应 深度 涉 
及 整个 观测 剖面 。 

对 降雨 量 和 日 JA .半月 月 尺度 各 层 土壤 含水 
量 进行 相关 分 析 ,发现 降 雨量 和 深度 较 浅 的 土壤 体 
积 含水 量 之 间 存 在 较 强 的 相关 关系 。 在 毛乌素 沙 
地 , 匡 丽 丽 等 "对 不 同时 间 尺 度 累积 的 降雨 量 和 沙 


层 渗 漏水 量 进行 相关 分 析 和 线性 拟 合 后 发 现 , 越 往 
深 处 渗 漏 水 量 对 降雨 的 响应 越 弱 。 廉 泓 林 等 ”以 
科尔沁 沙 地 南 缘 的 榜 子 松 和 柠 条 固沙 人 工 林 为 研 
究 对 象 ,进行 土壤 水 分 的 时 空 变化 特征 研究 ,得 到 
土壤 含水 量 对 大 气 降 雨 的 响应 随 着 土 层 深度 的 增 
加 而 减弱 的 结果 ,与 本 文 研究 结果 一 致 。 此 外 ,本 
试验 发 现在 一 定 范围 内 下 渗 深 度 随 降雨 量 的 增 大 
而 增加 , 当 降雨 量 大 于 50 mm 时 ,能 保证 对 2.0 m 深 
度 土壤 水 分 的 补给 。 
3.3 降水 再 分 配 过 程 模 拟 与 分 析 

Hydrus-1D 模 型 已 成 功 应 用 于 不 同类 型 林 分 的 
土壤 水 分 动态 模拟 中 ,包括 杨 树林 ,热带 季 雨 林 、 药 
漠 绿 洲 等 。 李 琦 等 3 利用 Hydrus-1D 模 拟 华北 平原 


Ez 
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小 麦 种 植 区 水 盐 运 移 ,结果 表明 浅 层 土壤 水 分 模拟 
误差 更 大 。 李 阳 明 ”在 南方 红壤 坡地 土壤 水 分 运 
移 模拟 中 也 得 到 模型 精度 在 空间 上 呈现 出 深层 高 
于 表层 的 特点 。 在 本 研究 8 个 观测 层 中 ,模型 对 0.4 
m 深 度 土壤 水 分 波动 模拟 较 其 他 深度 拟 合 情况 差 ， 
且 深 层 土壤 含水 率 的 模拟 精度 高 于 浅 层 土壤 。 浅 
层 土壤 水 分 运 移 模拟 因 受 外 界 因素 的 剧烈 影响 , 导 
致 比 深层 土壤 产生 更 大 的 误差 。 渗 透 .蒸腾 和 莹 发 
等 水 分 通 量 的 频繁 交换 可 能 造成 根系 主要 分 布 层 
产生 较 大 的 模拟 误差 “|。 

本 研究 应 用 参数 优化 的 Hydrus-1D 模型 模拟 了 
樟 子 松林 地 在 生长 季 的 水 分 动态 变化 ,结果 表明 ， 
整体 剖面 尺 为 0.876,8 个 观测 层 尺 值 范围 处 于 
0.61~0.85 ,相对 误差 绝对 值 和 均 方 根 误 差 值 的 范围 
分 别 为 0.0017~0.063 cm°- cm” F1 0.0061~0.0096 cm- 
cm ，。 李 冰冰 等 “指出 多 数 基于 Hydrus-1D 的 模 
拟 ,决定 系数 在 0.59~0.84, 均 方 根 误差 在 0.015~ 
0.063 cm em 。 本 研究 中 评价 指标 均 落 于 或 部 分 
优 于 上 述 范围 ,表明 模型 能 较 好 地 模拟 研究 区 土壤 
水 分 的 动态 变化 特征 ,模型 精度 较 高 。 此 校正 模型 
可 用 于 预测 降雨 后 土壤 水 分 的 分 布 与 运动 ,为 科 尔 
沁 沙 地 生态 修复 提供 科学 依据 。 


4 结论 

CD 种 植 樟 子 松 进行 生态 修复 后 改变 了 研究 区 
降雨 水 分 的 分 布 , 裸 沙 地 2 mm 深度 处 的 渗 漏 量 占 降 
雨量 的 44.16% , 樟 子 松林 地 深层 渗 漏 量 仅 占 降雨 量 
的 0.7%。 樟 子 松林 地 2 mm 处 深层 渗 漏 主要 发 生 在 
5 一 7 月 ,但 试验 期 间 林 地 仍 可 维持 4 水 转换 ,并 保证 
降雨 水 分 对 深层 渗 漏 或 者 地 下 水 的 补给 。 

(2) 监测 期 内 ,林地 0.4 m 深 度 以 下 土壤 水 分 对 
小 雨 无 响应 ,土壤 水 分 对 中 雨 的 响应 深度 可 达 1 m, 
对 大 雨 和 暴雨 的 响应 深度 涉及 整个 观测 剖面 ,1.2 m 
深度 土壤 体积 含水 量 在 全 年 均 高 于 其 他 深度 的 土 
壤 仿 水量。 随 着 土壤 深度 的 增加 ,林地 土壤 水 分 波 
动 幅度 呈现 减 小 的 趋势 ,0.4 m 深 度 土壤 水 分 受降 雨 
影响 强烈 。 

(3) 林地 降雨 量 和 深度 较 浅 的 土壤 体积 含水 量 
之 间 存 在 较 强 的 相关 关系 ;周期 为 周 .半月 的 累计 
降雨 量 与 各 层 土壤 体积 含水 量 显著 相关 , 当 周 期 为 
月 时 其 相关 性 整体 减弱 ,而 周期 为 日 时 其 相关 性 只 
在 较 浅 的 沙 层 中 体现 。 

(4) Hydrus-1D 模 型 可 准确 地 模拟 研究 区 土壤 


水 分 运动 , 且 深 层 模拟 精度 高 于 浅 层 。 
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Dynamic changes in soil moisture and its response to rainfall in Pinus sylvestris 
var. mongolica plantation in Horqin Sandy Land 


JIJI Jiamen", CHENG Yiben', CHEN Linglong, WAN Pengxiang, ZHANG Yihui, 
YANG Wenbin, BAI Xuying, WANG Tao 
(1. School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China; 2. Institute of 
Desertification Studies, Chinese Academy of Forestry, Beijing 100091, China; 3. Inner Mongolia Low Coverage Sand 
Control Technology Development Co. , Ltd. , Hohhot 010010, Inner Mongolia, China) 


Abstract: After ecological restoration in Horqin Sandy Land, the soil moisture dynamics of the forest and its 
response to rainfall affect whether Pinus sylvestris var. mongolica(PSM) can be continuously used for the 
ecological restoration in the same type of area. To study the characteristics of soil moisture dynamics after 
vegetation restoration on the southern edge of Horqin Sandy Land, this study used the methods of in situ 
observation and numerical simulation and corrected the HYDRUS-1D model based on the measured data of soil 
moisture to explore its relationship with rainfall. We obtained the following results: (1) The regional moisture 
distribution was changed by the PSM plantation, deep soil recharge at 2.0 m in bare sand accounted for 44.16% 
of the annual rainfall, whereas deep soil recharge at 2.0 m in PSM land accounted for only 0.7% of the annual 
rainfall. (2) Soil moisture below a depth of 0.4 m had no response to light rain, but the response depth of soil 
moisture to moderate rain could reach 1.0 m. The response depth to heavy rain and rainstorm involved the entire 
observation profile. With increasing soil depth, a decreasing trend was observed in the variation in moisture 
fluctuation. (3) There was a strong correlation between rainfall and volumetric soil water content at shallow 
depth, and the cumulative rainfall at weekly and semi- monthly intervals was significantly correlated with the 
volumetric soil water content of each layer(P«0.05). Rainfall >50 mm can ensure the supply of soil moisture 
within 2.0 m. (4) The determination coefficient of the model was between 0.61 and 0.85, and the root mean 


3 


square error ranged from 0.0061 to 0.0096 cm ?-cm ^. The accuracy of the deep layer simulation was higher than 
that of the shallow layer. These results have important implications for rain- fed vegetation afforestation in the 
Horqin Sandy Land. 

Keywords: rainfall infiltration; soil moisture; deep soil recharge; Pinus sylvestris var. mongolica; Horqin San- 


dy Land 


